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A high statistics measurement of the A+ lifetime from the Fermilab fixed-target FOCUS photopro- 
duction experiment is presented. We describe the analysis technique with particular attention to the 
determination of the systematic uncertainty. The measured value of 204.6±3.4 (stat.)±2.5 (syst.) fs 
from 8034 ± 122 A c — > pKn decays represents a significant improvement over the present world av- 
erage. 

PACS numbers: 13.30.Eg, 14.20.Lq, 14.65.Dw 



Experimental measurements of charm particle lifetimes 
have been used in the study of strong interaction physics. 
The measurements provide some guidance for theoret- 
ical calculations of non-perturbative strong interaction 
processes. The steady improvement in the precision of 
the measurements has not only helped to improve our 
theoretical understanding of strong interactions, but also 
to help stimulate the development of better theoretical 
tools. These have progressed from the spectator model 
to various quarks models and currently to Heavy Quark 
Expansion methods . These calculational tools are the 
same or similar to those used in other areas, for example 
to determine the size of the V u b CKM element through 
inclusive semileptonic B decays [0. More precise mea- 
surements of all of the charm particle lifetimes will help 
continue this process of improvement and extension of 



applicability. 

Precise charm lifetime measurements are now begin- 
ning to emerge from e + e~ collider experiments Q. 
The effects of lifetime and vertex resolution are also im- 
portant in mixing and CP violation measurements [j5], [| . 
It is crucial to have accurate lifetime measurements from 
fixed-target experiments to act as a standard to evalu- 
ate any relative systematic differences. The A+ lifetime 
presented in this paper represents the most accurate mea- 
surement of this quantity to date and is a significant im- 
provement over the present world average. 

The data used were collected by the FOCUS collab- 
oration in the 1997 fixed-target run at Fermi National 
Accelerator Laboratory. The FOCUS spectrometer is an 
upgrade of the spectrometer used in the E687 photopro- 
duction experiment M. The vertex region consists of 
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FIG. 1: pKn invariant mass plot for data (points) fitted 
with a Gaussian signal and quadratic background (solid line) . 
The shaded area indicates the fitted level of background. The 
vertical dotted lines indicate the signal and sideband regions 
(see text) used in the lifetime analysis. 



four BeO targets and 16 planes of silicon strip detec- 
tors (SSD). Two of the SSD planes were placed imme- 
diately downstream of the second target, and two im- 
mediately downstream of the fourth (most downstream) 
target. Momentum analysis was made possible by the use 
of 5 multiwire proportional chambers and two magnets 
with opposite polarities. Hadronic particle identification 
was achieved using three multicell threshold Cerenkov 
counters The data for this measurement were taken 
using a photon beam with average energy of ~ 180 GeV 
for triggered events. 

The A+ — > pK ~7r + Jl3[ candidates are reconstructed 
using a candidate driven algorithm which is highly effi- 
cient for all decays including short lived ones. All pK~n + 
candidates are tested to see if they form a vertex with a 
confidence level greater than 1%. The candidate A+ mo- 
mentum vector is then projected to search for a produc- 
tion vertex with one or more tracks. As many tracks as 
possible are included in the production vertex so long as 
the vertex confidence level is larger than 1%. The pro- 
duction vertex is required to be within one of the four 
targets. The separation L between the production and 
decay vertices is required to be larger than 6a i where ol 
is the error on L calculated on a candidate-by-candidatc 
basis. In addition, each track in the pK~Tr + candidate 
combination must also satisfy the appropriate Cerenkov 
particle identification criteria. 

The pKn invariant mass plot for data is shown in 
Fig. g. The fit shown uses a Gaussian signal and a 
quadratic background function which yields 8034 ± 122 
reconstructed A c decays. The lifetime analysis uses pK tt 
candidates within the signal and symmetric sideband re- 
gions as shown in the figure. All three regions are 4cr m 
wide and the centers of the sideband regions are lo- 
cated ±6(T m from the mean of the fitted Gaussian, where 
a m — 8.2 MeV/c 2 is the width of the fitted Gaussian. 

For the lifetime analysis we use the reduced proper 
time, t' = (L — 6ctl)//?7C where /3j = pA c /rriA c and 
require it to be less than 1 ps to reduce long-lived back- 
grounds. This requirement was already made for the data 
shown in Fig. |l|. The use of the reduced proper time en- 



sures that only a small acceptance correction to the life- 
time distribution is needed. The average proper time res- 
olution for this decay sample (42 fs) is small enough com- 
pared to the lifetime to use a binned likelihood method 

1- 

The t' distributions for the decays in the signal and 
sideband regions are binned into two separate histograms 
from 0-1 ps in 20 fs bins. The observed number of de- 
cays in the 4 th t' bin is s, for the signal region and 6j for 
the sideband region. The t' distribution of the sideband 
region is used as a measure of the lifetime distribution 
of background events in the signal region. Thus the ex- 
pected number of decays in the « th t' bin of the signal 
region is given by: 
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The likelihood that is maximized in the fit is given by 

(aB) Nb e~ aB 



(1) 



Likelihood = 



n ■ e 



(2) 



where S is the total number of signal events and B is the 
total number of background events in the signal region 
and S + B = Esj. The total number of events in the 
sideband region is Nf, — and a is the ratio of the 
number of events in the sideband region to the number 
of background events in the signal region. The value of a 
is obtained from the fit to the invariant mass distribution 
and is very close to 2. B and r are the fit parameters. 

The effects of geometrical acceptance, detector and re- 
construction efficiencies, and absorption are given by the 
f(t') correction function. The f(t') is determined us- 
ing a detailed Monte Carlo (MC) simulation of the ex- 
periment where the production (using Pythia jlC[]) was 
tuned so that the production distributions for data and 
MC matched. Note that only the shape of the f(t') func- 
tion is important and it is obtained by dividing the ob- 
served MC t' distribution by a pure exponential with the 
MC generated lifetime. The f(t') distribution is shown 
in Fig. |(a). 

Using the likelihood function given above we obtained 
a fitted lifetime of 204.6 ±3.4 fs. The lifetime distribution 
of all decays in the signal region is shown in Fig. ||(b) to- 
gether with the fit and the level of background contained 
in the signal region. 

Detailed studies were performed to determine the sys- 
tematic uncertainty in this measurement. 

The uncertainty in the absolute time scale was inves- 
tigated by studying the absolute length and momentum 
scales in the experiment. For the length scale, compar- 
isons were made between measurements of the distances 
between silicon planes in the target region. The values 
obtained using vertex positions in the data with the stan- 
dard vertexing code agree well with those obtained us- 
ing precision instruments. The absolute momentum and 
mass scales were checked by comparing the reconstructed 
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FIG. 2: (a) The f(t') correction function. Deviation from 
a fiat line indicates the correction from a pure exponential; 
(b) The lifetime distribution for all decays in the data sig- 
nal region (points), and the fit (histogram). The shaded dis- 
tribution shows the lifetime distribution of the background 
component in the signal region; (c) The lifetime distribution 
for A c decays (points), i.e. the sideband subtracted and f(t') 
corrected yield. The line is a pure exponential with the fit- 
ted lifetime and the shaded region gives the background. An 
arbitary yield scale is used because of the particular normal- 
ization of fit'). 



masses of charm and strange mesons and hyperons with 
established values. Our studies showed no evidence of 
any scale offset, but due to the limited statistical pre- 
cision of these comparisons we assign an uncertainty of 
±0.11% to the absolute time scale. 

The backgrounds are composed of a non-charm and 
a charm component; these two background components 
are approximately equal in our sample and fairly evenly 
distributed across the signal and sideband mass regions. 
The level and lifetime distribution of the background in 
the signal mass region is assumed to be well represented 
by symmetric mass sidebands close to the signal region. 
The uncertainties that arise because of these assumptions 
were determined by a large number of studies. 

The contamination from D + — > K^ir + n + , D + — > 
K~K + n + and Df — > K~K + n + decays misidentified as 



decays were determined in our sample. We loos- 



ened the Cerenkov requirements on the data and used the 
MC efficiencies to extrapolate to tighter particle identi- 
fication criteria. From this we found the above three 
decays respectively contribute 0.5%, 1.3% and 2.7% of 
the total background in the signal region. The small 
contribution of these reflection backgrounds and the fact 
that they are distributed fairly uniformly across the sig- 
nal and sideband mass regions mean they give rise to 
insignificant uncertainties. This was verified in a test by 
explicitly eliminating them by cutting out the appropri- 
ate mass regions. Using variations in particle identifica- 
tion and vertexing selection to significantly change the 
signal/background ratio also showed no significant un- 
certainties. 

The background lifetime uncertainty was further inves- 
tigated by using symmetric sidebands of different widths 
(4-16<r m ), and located at different separations from the 
signal region (±4 to ±16<r m ). The effect of using only the 
low or only the high mass sideband was also studied. The 
effect of having the fit parameter B truly free by eliminat- 
ing the background term in the likelihood (second term 
in Eq. (|^)) was studied and found to be inconsequential. 
Note that the results of the pKir mass fit are only used 
in the background term in the likelihood. 

Finally, an independent analysis which did not rely on 
knowledge of the background lifetime distribution was 
performed. In this analysis the data were split into 
twenty 50 fs wide reduced proper time bins from 0-1 ps. 
The number of A+ — > pK~n + decays in each bin was 
determined in a mass fit and the yields fitted to an expo- 
nential decay distribution modified by a f(t') correction 
function. This fit') function was obtained separately 
for this analysis from the MC, doing the same split into 
twenty time bins and fitting the mass distributions for 
each MC bin. This f(t') correction function agrees well 
with that obtained in the standard analysis method. 

From these studies we assign a background systematic 
uncertainty of ±0.77%. 

Uncertainties in the f(t') correction include uncertain- 
ties from the geometrical acceptance, the detector and 
reconstruction efficiencies, the production model, the ab- 
sorption cross-sections, and the decay dynamics. 

With our chosen selection criteria, the fit') correction 
reduces the fitted lifetime by 1.19%. A number of studies 
were performed to study the uncertainty in this correc- 
tion. Since the correction function is obtained from MC 
simulations, care was taken to ensure that this simulation 
correctly reproduces a very large number of data distri- 
butions. In particular the MC reproduces the data A+ 
longitudinal and transverse momenta, the multiplicity of 
the production vertex, and the decay length and proper 
time resolutions. A sensitive check of the acceptance and 
efficiency part of the MC correction was done using high 
statistics K s — ► tt + it~ decays. Short-lived Kg decays 
were reconstructed using the same analysis methods in 
the same decay region as the A+ decays. Since the K s 
lifetime is well known we can determine the f(t') correc- 
tion in data and compare it to that obtained in our MC 
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simulation. The agreement is excellent but was limited 
by both data and MC statistics to a sensitivity of ±2% of 
the correction. Using this as the level of the uncertainty 
in the fit') correction, we can assign a systematic uncer- 
tainty due to this correction of ±0.83%. Possible time 
dependent systematic effects were looked for by splitting 
the data into different time periods and comparing the 
fitted lifetimes. We also compared the separate fitted 
lifetimes for decays originating from each of the four tar- 
gets. No systematic uncertainties were found in these 
two comparisons. 

Our limited knowledge of the production and decay 
of the A+ could contribute to a systematic uncertainty. 
This was studied using different MC simulations where 
the production parameters and the resonance substruc- 
ture of the decay were varied over reasonable ranges. 
Production systematics were also studied by splitting the 
data into different bins of longitudinal and transverse 
momenta, primary vertex multiplicity, and by comparing 
the fitted lifetimes for particles and anti-particles. We 
assign a systematic uncertainty of ±0.38% due to our 
limited knowledge of A+ production and decay. 

In order to use the reduced proper time we must be 
able to correctly model our proper time resolution. This 
was verified by comparing the distributions for data and 
MC and by studying splits of the data sample that can 
be sensitive to resolution effects. The data were split into 
bins of proper time resolution and reconstructed invari- 
ant mass. Variations of the proper time bin width from 
10 to 100 fs were also studied as was changing the fitted 
range from 0-0.6 ps to 0-1.4 ps, and from 0-1 ps to 0.2- 
1 ps. We assign a systematic uncertainty of ±0.12% to 
the lifetime due to resolution uncertainties. 

The systematic uncertainty due to absorption of the 
A+ and daughter particles was studied by varying the 
charm interaction cross-section by 100% and the daugh- 
ter particle interaction cross-sections by 50% in the MC. 
It was also studied by comparing the lifetimes of decays 
occuring inside and outside of the target, and by compar- 
ing the lifetimes for decays where the A+ was produced in 
the upstream half of each target with those produced in 
the downstream half of the same target. We determined 
a systematic uncertainty of ±0.23% due to absorption. 

Contributions to the systematic uncertainty are sum- 
marized in Table [j]. Taking contributions to be uncorre- 
cted we obtain a total systematic uncertainty of ±1.23% 
or ±2.5 fs. 



We have measured the A+ lifetime to be 204.6 ± 
3.4 (stat.) ± 2.5 (syst.) fs using 8034 ± 122 A c -> pKn 
decays from the Fermilab FOCUS photoproduction ex- 
periment. This measurement represents a significant im- 
provement in accuracy and special care was taken to in- 
vestigate and properly quantify possible systematic un- 
certainties. Table [ij compares our measurement with 
previous recent published results. The difference between 
this measurement and the measurement from the CLEO 
e + e~ experiment may point to the emergence of possible 
relative systematic effects [Il2l . Any such systematic dif- 



TABLE I: Contributions to the systematic uncertainty. 



Contribution 




Systematic (%) 


Time scale 




±0.11 


Backgrounds 




±0.77 


Acceptance 




±0.83 


Production 




±0.38 


Resolutions 




±0.12 


Absorption 




±0.23 


Total 




±1.23 


TABLE 11: Comparison 


of recent A+ lifetime measurements. 


Experiment 


Type 


r(A+) fs 


E687 f] 


FT 


215 ±16 ±8 


SELEX |ll[ 


FT 


198.1 ±7.0 ±5.6 


CLEO 11.5 § 


e + e~ 


179.6 ±6.9 ±4.4 


FOCUS (this result) 


FT 


204.6 ±3.4 ±2.5 



ference would be important to resolve given the number 
of recent and future mixing and CP-violation measure- 
ments that rely on accurate knowledge of lifetime distri- 
butions. 

We wish to acknowledge the assistance of the staffs 
of Fermi National Accelerator Laboratory, the INFN of 
Italy, and the physics departments of the collaborating 
institutions. This research was supported in part by the 
U. S. National Science Foundation, the U. S. Department 
of Energy, the Italian Istituto Nazionale di Fisica Nucle- 
are and Ministero delPUniversita e della Ricerca Scien- 
tifica e Tecnologica, the Brazilian Conselho Nacional de 
Desenvolvimento Cientffico e Tecnologico, CONACyT- 
Mexico, the Korean Ministry of Education, and the Ko- 
rean Science and Engineering Foundation. 



[1] G. Bellini, I. I. Y. Bigi, and P. J. Dornan, Phys. Rept. 

289, 1 (1997). [8J 



[2] |M. B. Voloshin, Phys. Lett. B515, 74 (2001). 

[3] G. Bonvicini et al., Phys. Rev. Lett. 82, 4586 (1999). 
[4] A. H. Mahmood et al., Phys. Rev. Lett. 86, 2232 (2001). 
[5] R. Godang et al., Phys. Rev. Lett. 84, 5038 (2000). 
[6] K. Abe et al., hep-ex/0111026. 

[7] P. L. Frabetti et al., Nucl. Instrum. Methods Phys. Res., 



ex/0108011 



lep- 



Sect. A 320, 519 (1992). 

J. M. Link et al., FERMILAB-Pub-01/243-E, 

to be published in Nucl. Instrum. Methods 
Sect. A. 

Phys. Rev. Lett. 70, 1755 (1993). 
Comput. Phys. Commun. 135, 238 



Phys. Res 
[9] P. L. Frabetti et al 
[10] T. Sjostrand et al. 

(2001). 
[11] A. Kushnirenko et al 



Phys. Rev. Lett. 86, 5243 (2001). 



5 



[12] H. W. K. Cheung. FERMILAB-Conf-01/351-E, ^ep- 



ex/0111050 



[13] The charge conjugate mode is implicitly implied unless 

otherwise stated. 
[14] The reduced proper time has the clock started at the 



minimum allowed time for each decay candidate. In the 
absence of other corrections, it does not matter where 
along the decay exponential one starts the clock, hence 
the reduced proper time also follows a pure exponential. 



